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Abstract 

The cosmic-ray energy spectrum above 10 18,5 eV is reported using the updated 
data set of the Akeno Giant Air Shower Array (AGASA) from February 1990 
to October 1997. The energy spectrum extends beyond 10 20 eF and the energy 
gap between the highest energy event and the others is being filled up with 
recently observed events. The spectral shape suggests the absence of the 2.7 K 
cutoff in the energy spectrum or a possible presence of a new component 
beyond the 2.7 K cutoff. 

PACS numbers: 98.70.Sa, 96.40.Pq, 96.40.De 
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How high the maximum energy of cosmic rays reaches is one of the most important 
problems in cosmic ray research. Detections of cosmic rays with energies above 10 20 eV 
]1]J2|] have given rise to much discussion regarding their origin. Many models have been 
proposed as source candidates of such high energy cosmic rays: active astrophysical objects 
||, decay products of much higher energy particles such as superheavy relic particles [|J or 
topological defects [|J, or cosmological gamma-ray bursts @ (see Ref. |7j for a recent review). 
If such high energy cosmic rays come from far outside our Galaxy, they interact with cosmic 
microwave background photons and cannot travel cosmological distances. This interaction 
causes a cutoff in the energy spectrum near 5 x 10 19 eV which is referred to as the Greisen- 
Zatsepin-Kuz'min (GZK) cutoff ||. Furthermore, the cosmic rays which have interacted 
form a "bump" just below the GZK cutoff energy |]9|-|ll|. The change in the spectral slope 



around 10 19 eV ("ankle") may arise from a transition from galactic to extragalactic sources. 
The investigation of these features in the energy spectrum is one of the most important 
scientific challenges. 

There are two techniques for detecting extensive air showers (EAS): widely spread surface 
arrays and atmospheric fluorescence detectors. Using these techniques, the energy spectrum 
of extremely high energy cosmic rays has been measured by many groups such as Volcano 
Ranch @, Haverah Park |3J, Sugar [TJ, Yakutsk [T§, Fly's Eye [0, and Akeno |I7|g3 



(only the Fly's Eye group has adopted the atmospheric fluorescence detector). While the 
energy spectrum obtained from these experiments coincide within ±15% in energy below 
~ 10 19 eV, the details of energy spectrum in the highest energy range is still inconclusive, 
mainly because of low statistics of their observed events. In this letter, we present the energy 
spectrum above 10 18 ' 5 eV A obtained from the Akeno Giant Air Shower Array (AGASA) f[8,19 



which currently has the largest exposure of any extremely high energy cosmic ray detectors. 

The AGASA array is the largest operating surface array, covering an area of about 
100 km 2 and consisting of 111 surface detectors of 2.2 m 2 area. Each surface detector is 
placed with a nearest-neighbor separation of about 1 km and the detectors are sequentially 
connected with pairs of optical fibers. All the detectors are controlled at detector sites 
through rapid communication with a central computer. The data acquisition system of 
AGASA was improved in December 1995 ||19|| . In a widely spread surface array like AGASA, 
the local density of charged shower particles at a specific distance from the shower axis is 
well established as an energy estimator |[2Q|| , since this depends weakly on variation in the 
interaction model, fluctuation in shower development and the primary mass. In the AGASA 
experiment, we adopt local density 5(600) at 600 m which is determined from fitting the 



lateral distribution of observed particle densities to an empirical formula |21| . This empirical 
formula is found to be valid for EAS with energies up to 10 20 eV |2"2] , |2"3| ] . The conversion 
relation from 5(600) to the primary energy is evaluated through the Monte Carlo simulation 
23 up to 10 19 eV by 



E — 2.03 x 10 So (600) eV, 

where So (600) is the S(600) value in units of m~ 2 for a vertically incident shower. Since an 
inclined air shower traverses more atmospheric depth than a vertical shower, S#(600) ob- 
served with zenith angle 9 must be transformed into So (600) at the vertical. This attenuation 
curve of S(600) has been formulated by Yoshida et al 
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The accuracy of event reconstruction has been evaluated through the analysis of a large 
number of artificial air shower events. These artificial events were simulated over a larger 
area than the AGASA area with directions sampled from an isotropic distribution. In 
this air shower simulation, the fluctuation on the longitudinal development of air showers, 
the resolution of the scintillation detectors, and statistical fluctuation of observed shower 
particles at each surface detector were taken into account. Only events with zenith angles 
smaller than 45° and with core locations inside the array area are used in the following 
analysis. Fig. [I] shows the fluctuation of energy determination for 10 19,5 e\^ (left) and 10 20 eV 
(right) showers with zenith angles less than 45°. The primary energy is determined with an 
accuracy of about ±30% and the proportion of events with a 50%-or-more overestimation 
in energy is about 2.4%. 

Energy uncertainty also arises from the following systematic errors. The first is uncer- 
tainty in measuring the particle density incident upon each detector. The number of incident 
particles is determined from the time width of a pulse, which is generated by decaying an 
anode signal of a photomultiplier tube exponentially with a time constant of about 10/xs 
and discriminated at a certain level (see ]H| for the details of the AGASA instruments). 
The variation in the amplifier gain and the decay constant are monitored in every run for 
detector calibration and their seasonal variations are within 2%. The second is uncertainty 
in the empirical formula of the lateral distribution function and in the attenuation curve of 
5(600). The energy uncertainty due to the limited accuracy on both of these is estimated 
to be ±20%, even if both factors shift the estimated energy in the same direction The 
third is uncertainty in the conversion formula of 5(600) into primary energy. Although 
this formula is not sensitive to interaction models or primary composition in each of sim- 
ulation codes HU, the systematic errors due to the differences in simulation codes are not 
quantitatively clear. 

In order to evaluate the systematic errors experimentally, we compare the AGASA spec- 
trum derived below with the Akeno spectrum which was accurately determined between 
10 14 ' 5 eV^ and 10 19 eV using the arrays with different detector spacing ]17|]. The Akeno spec- 
trum fits very well with extrapolation of those obtained from direct measurement on balloons 
and satellites, and with the Tibet result |26[ obtained through the observation of the shower 
at the height of its maximum development. The difference between the present AGASA 
and Akeno spectra is about 10% in energy at 10 l8 5 eV. In addition, the difference among 
spectra obtained from the Fly's Eye, Yakutsk, Haverah Park, and AGASA experiments is 
within 30% in energy in spite of quite different methods for determining the primary energy. 
Therefore, the total systematic error in the AGASA energy estimation is estimated to be 
within 30%, and the primary energy of the highest energy event of AGASA, for example, is 
estimated to be in the range (1.7 - 2.0) xl0 20 e^. 

The effective area of AGASA has been calculated from the simulation of artificial air 
shower events. The energy spectrum in this simulation was assumed to be E~ 3 , and the 
reconstruction uncertainty in energy estimation was also taken into account. Although the 
effective area depends weakly on the spectral index, this dependence is negligible when 
compared with other ambiguities like energy resolution. The total exposure of AGASA is 
obtained by multiplying the effective area and the observation time of each branch for each 
epoch. Above 10 19 eV, this exposure is constant and is 2.6 x 10 16 m 2 sr s, which is about five 
times as large as that in our previous paper ]23| (cf. ~ 0.5 x 10 16 m 2 sr s of the stereo Fly's 
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Eye exposure |16j and ~ 0.7 x 10 m sr s of the Haverah Park exposure ||13|| ). However, 



the exposure below 10 18 5 eV depends strongly on the primary energy. Since this energy 
dependence causes systematic errors in the energy spectrum derivation, only events with 
energies above 10 18,5 eV^ are used for the energy spectrum in this letter. From February 1990 
to October 1997, 3847, 461 and 6 events were observed with energies above lO 18 ' 5 ^, 10 19 eV^ 
and 10 20 eV respectively. 

The energy spectrum observed with AGASA is shown in Fig. |2|, multiplied by E 3 in order 
to emphasize details of the steeply falling spectrum. Error bars represent the Poisson upper 
and lower limits at 68% and arrows are 90% CL. upper limits. Numbers attached to points 
show the number of events in each energy bin. The dashed curve represents the spectrum 
expected for extragalactic sources distributed uniformly in the Universe, taking account of 



the energy determination error flTT 



First, we examine whether the observed energy spectrum could be represented by a single 
power law spectrum (oc E^ 11 ). The optimum spectral index 71 is derived from the maximum 
likelihood procedure comparing the observed and expected number of events in each energy 
bin. This procedure is same as described in Yoshida et al. |[25|| . The maximum likelihood 
procedure for a single power law spectrum results in 71 = S-OS^oj 8 ; the likelihood significance 
of 71 is only 0.051. If only events with energies below 10 19 eV are considered, 71 (E < 
10 19 e\^) = 3. 23to;}2 is obtained which is consistent with the spectral index, 3.16 ± 0.08, 



determined from the Akeno experiment [17 



Next, a broken energy spectrum is examined with the same procedure, 
energy spectrum is assumed to be 

dJ__{n {E/E a y 10 10 18 - 5 eV < E < E a 
dE~ \K(E/E a )^ 2 E a < E 



The broken 



where 70 and 72 are indexes below and above a bending (ankle) energy E a , and 70 is fixed 
to be ji(E < 10 19 e^) = 3.16 determined from the Akeno experiment [I7|. The most 
probable parameters are obtained at E a = 10 19 ' 01 e\^ and 72 = 2.78ioJ3, where the likelihood 
significance is found to be 0.903. This is also consistent with the results of 2.8±0.3 at energies 
above 10 18 ' 8 eV determined from the Akeno experiment [17] and of 2.3 
previous paper 



l+0 i above 10 



19.0 



in the 
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Furthermore, the energy spectrum presented here extends up to higher energies than the 
previous results |L7,2S[]; six events were observed above 10 20 eV. If the real energy spectrum 
is that shown in Fig. |2] as the dashed curve, the expected number of events above 10 20 eV 
is less than one, taking account of the energy resolution. The energy spectrum is therefore 
more likely to extend beyond 10 20 eV without the GZK cutoff. However, it is also worth 
noting that the observed energy spectrum suggests a small deficit just below 10 20 e\^, whose 
significance is not compelling because of the uncertainty in 72 estimation. This deficit may 
imply another component above the GZK cutoff energy. In either case, sources of the most 
energetic cosmic rays must be located within a few tens of Mpc from our Galaxy |TT . 
The arrival directions of six 10 20 eV events are shown in Fig. [5| Within the accuracy of 
arrival direction determination (1.6° above 4 x 10 19 eV), no 10 20 eV events coincide with 
possible candidates from the second EGRET sources |27j or the extragalactic radio sources 
with redshift z < 0.02 p8| . Our previous result for cosmic-ray arrival directions has been 



reported in Hayashida et al. p| and the new results are under preparation. 
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The fact that the energy spectrum extends beyond 10 20 eV and no 10 20 eV events coincide 
with nearby active astrophysical objects leads highest energy cosmic-ray physics into a much 
more exciting stage. The next generation experiments such as the Telescope Array pO| , |3T 



High Resolution Fly's Eye |32],[33|, and Auger [@,^] projects will solve the puzzle of the 
highest energy cosmic rays. 

In conclusion, the cosmic-ray energy spectrum extends beyond 10 20 eV. No candidate 
sources are found in the directions of six 10 20 eV events, while their sources must be closer 
than 50 Mpc. The possible deficit around 10 20 eV is a notable area in which to search for 
origin of the highest energy cosmic rays. Detailed discussion with the AGASA data will be 
published elsewhere. 

We are grateful to Akeno-mura, Nirasaki-shi, Sudama-cho, Nagasaka-cho, Ohizumi-mura, 
Tokyo Electric Power Co. and Nihon Telegram and Telephone Co. for their kind cooperation. 
The authors are indebted to other members of the Akeno group in the maintenance of the 
AGASA array. The authors thank Jamie Holder and Michael Roberts for their valuable 
suggestion on the preparation of the manuscript. 
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FIGURES 



FIG. 1. Fluctuation of energy determination for 10 195 eU (left) and 10 20 eU (right) showers 
with zenith angles less than 45°. 

FIG. 2. Energy spectrum observed with AGASA. The vertical axis is multiplied by E 3 . Error 
bars represent the Poisson upper and lower limits at 68% and arrows are 90% C.L. upper limits. 
Numbers attached to points show the number of events in each energy bin. The dashed curve 
represents the spectrum expected for extragalactic sources distributed uniformly in the Universe, 
taking account of the energy determination error [11]. 

FIG. 3. Arrival directions of six 10 20 eU events on the Galactic coordinates. The shaded regions 
indicate the non-observable celestial regions due to the zenith angle cut of < 45°. The equatorial 
and supergalactic planes are also shown. 
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FIG. 1. Fluctuation of energy determination for 10 19 - 5 eF (left) and 10 20 eV (right) 
with zenith angles less than 45°. 
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FIG. 2. Energy spectrum observed with AGASA. The vertical axis is multiplied by . Error 
bars represent the Poisson upper and lower limits at 68% and arrows are 90% C.L. upper limits. 
Numbers attached to points show the number of events in each energy bin. The dashed curve 
represents the spectrum expected for extragalactic sources distributed uniformly in the Universe, 
taking account of the energy determination error [11]. 
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FIG. 3. Arrival directions of six 10 eV events on the Galactic coordinates. The shaded regions 
indicate the non-observable celestial regions due to the zenith angle cut of < 45°. The equatorial 
and supergalactic planes are also shown. 
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